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Synchronization of an array of spin torque nano oscillators in pe- 
riodic applied external magnetic field 
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Abstract - Considering an array of spin torque transfer nano oscillators (STNOs), we have 
investigated the synchronization property of the system under the action of a common periodically 
driven applied external magnetic field by numerically analyzing the underlying system of Landau- 
Lifshitz-Gilbert-Slonczewski (LLCS) equations for the macro-spin variables. We find the novel 
result that the applied external magnetic field can act as a medium to induce synchronization of 
periodic oscillations, both in-phase and anti-phase, even without coupling through spin current, 
thereby leading to the exciting possibility of enhancement of microwave power in a straightforward 
way. 



INTRODUCTION. - It has been increasingly re- 
alized that the STNO can be developed into a useful de- 
vice for submicron microwave generation due to its varied 
oscillatory properties as a function of injected spin cur- 
rent and appUed magnetic fields [T1|31[T3].A rich variety of 
dynamics has been found in a single STNO by applying 
a periodically varying spin current [11] or applied mag- 
netic field [5Tj. Even though an individual STNO, which 
is an interesting nonlinear dynamical system, has prac- 
tical limitations due to its low output oscillation power 
[1] , nevertheless it has the interesting possibility of gener- 
ating increased power though the process of synchroniza- 
tion with other similar STNOs under appropriate coupling 
|15| . Synchronization of nonlinear oscillators both in the 
periodic and chaotic regimes are topics of intense discus- 
sion in recent times in the nonlinear dynamics literature 
[9l[T0]. For STNOs it has been shown that mutual phase 
locking phenomenon can be realized leading to increased 
output power through synchronization in both series and 
parallel architectures by electrical coupling of spin cur- 
rent with and without delay [T51 [H] .Although there are 
various ways of achieving synchronization, even four hun- 
dred years ago Huygens observed the same by applying a 
common load to a system of two pendulum clocks. Re- 
cently Nakada, Yakata and Kimura [16] have also pointed 
out the interesting possibility of noise induced synchro- 
nization of a pair of STNOs under noisy current injection 



for controlling the output power in the array of STNOs. 
This particular study shows that an application of com- 
mon Gaussian white noise to a system of two uncoupled 
STNOs acts as a medium for them to induce synchroniza- 
tion. 

In this paper, we present an alternative straight forward 
approach of synchronization of the macro-spin dynamics 
of an array of STNOs through the simple mechanism of 
applying a common driven external periodic applied mag- 
netic field in the presence of a dc magnetic field and dc 
spin current, without any further coupling. The applied 
periodic magnetic field acts as a medium through which an 
effective coupling occurs between the macro-spins of the 
free layers of the different STNOs leading to an effective 
synchronization involving either in-phase locked or anti- 
phase locked dynamics. We also present a phase diagram 
for a system of two STNOs delineating the in-phase and 
antiphase synchronization and desynchronization regions 
in the dc magnetic field vs dc spin current plane. We also 
present the results for 4 and 100 STNOs to confirm the 
possibility of synchronization of a large number of STNOs, 
even in the presence of mismatches in system parameters. 

Nonlinear dynamics of a Spin Torque Nano Os- 
cillator. — The dynamics of the magnetization in the 
free layer [Sl[7] of an STNO is known to be well described 
by the LLCS equation for the normalized spin vector S, 
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which is a classical unit vector, S = Sxi + Syj + Szk, 
I^P = 1. It is given by the nonlinear evolution equation, 
see for example ref. [5], 



We assume the spin polarization vector Sp = (1,0,0) 
and that the applied external magnetic field consists of dc 
and ac components. Therefore, 



dS_ 
It ' 

where 
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Heff ='Hi + nSxi — 47rS'oS'2fc, 71(1) = hdc + hac cos{rtt). 

(3) 
Under a stercographic projection, eq. ([T]) can be equiv- 
alently rewritten using the complex scalar function (18j 
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In the above a is the density of the spin current, defined 
as a = hT]j /2SoVe, ry is the area of cross section of the 
pillar, j is the charge current density, 5*0 is the saturation 
magnetization, V is the volume of the fixed layer. It is 
to be noted that the the density of the spin current is 
proportional to the density of the charge current injected 
into the fixed layer. Note that with the above definition 
a can be expressed in the unit of Oersted [TT] . Hence by 
varying the charge current one can vary the spin current 
injected into the free layer. 
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Fig. 1: The largest Lyapunov exponent (A) is plotted in the a- 
hdc plane for a single STNO for an oscillating external magnetic 
field of strength hac = 10 Oe and frequency £7 = 15 GHz with 
anisotropy parameter k = 45.0 Oe along the in-plane axis. 
The colour coding clearly distinguishes periodic and chaotic 
regimes. 

Considering homogeneous spin states on the ferromag- 
netic free layer, we can treat H^xchange to be zero. We 
assume an easy plane anisotropy, Hanisotropy = {kSx, 0, 0) 
where n is the strength of the anisotropy. For a nano pil- 
lar geometry, we assume H^^mag = — 47rS'o(0, 0, 5^) and 
Happiied = (^(i),0, 0) so that the external magnetic field 
is applied along the x-direction. In our following analysis, 
as illustration, we choose the parameters typically as that 
of a permalloy thin film. So we have fixed the Gilbert 
damping parameter as A = 0.02, saturation magnetiza- 
tion as 47rS'o = 8.4 kOe, and the gyromagnetic ratio as 
7 = 0.01767 Oc-i ns-i. 



^(i) 



_ CJ + CJ 

^" l + |a;P ' " 



"--'w 



bx + IJy 
1 + 5. ' 
.{lJ — id) 
"*l+|wP 



(4) 



5z = 



1-lc^l 



7(0 



f)(l--^) 



UJ + UJ* 



2 ' 1 



- OJU! 



.(1-..^) 



z747rS'o 



1 



1 



(5) 



One may note that the effect of applied magnetic field 
and the spin current terms occur together in the first term 
in the right hand side of eq. ([Sj . So one may expect the ef- 
fect of spin current and magnetic field to complement each 
other. It has been shown [THIIH] that in the case of ho- 
mogeneous magnetization in the free layer, where the spin 
S" is a function of time only, S = S(t), LLCS equation 
can exhibit interesting bifurcation and chaos scenario de- 
pending on the strength and form of the applied external 
magnetic field and injected current. In fig. [U we present a 
two parameter (a— hj,„) phase- like diagram characterized 
by the value of the largest Lyapunov exponent which es- 
pecially distinguishes periodic and chaotic regimes. This 
is obtained numerically by integrating eq. (P) or eq. ([5]) of 
a single STNO using embedded variable step size Runge- 
Kutta method in the current (a) Vs applied magnetic field 
(hdc) plane with anisotropy strength k = 45 Oe, where the 
other parameters are chosen as given in the figure caption. 
Note that the spin current (a) and the applied magnetic 
field (hdc) are in units of Oersted(Oe). 

All our numerical analysis have been done on a 32-node 
cluster facility with 8 digit accuracy in the normalization 
of the magnetization vector S. In fig. [1] the yellow/green 
regions are periodic, where all the Lyapunov exponents 
are zero or negative. In the red region one of the Lya- 
punov exponents (among the three) is greater than zero, 
that is positive as their values are encoded in color palette. 
One may note that in the lower ranges of magnetic field 
and current strength limit cycles/periodic oscillations arc 
present, whereas in the higher strength regions the spins 
are in chaotic motion. One can also check that typical 
transition from periodic to chaotic states occur through 
well defined routes such as period doubling bifurcations 

m. 
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Synchronization of STNOs in periodic applied magnetic field 



Synchronization of coupled STNOs. The 

schematic representation of our model for an array of two 
STNOs (see fig. [2]) in the presence of a common applied 
magnetic field 



Happ = [hdc + hac cos Vtt, 0, 0), 



(6) 



where each one of the STNOs is having a separate dc 
charge current source and the synchronized output cur- 
rent from each STNO as taken in its free layer is added 
as shown in fig. [51 Considering a system of two electri- 
cally uncoupled STNOs placed in a common magnetic 
field®, the magnetization dynamics is governed by the 
LLCS equations for the magnetization vectors 81,82, 
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Fig. 2: Schematic representation of an array of two STNOs 
placed in the oscillatory external magnetic field (Happ)- The 
-|- sign stands for summing up of two output currents. 



-T^ = -'ySi X i?igfj + XSi X — i - jaSi X (Si X Sp), 



dt 

^ = -7& X J?2„ff + A& X ^ - jaS2 X (& X Sp) (7) 
at '^^ 



dt 

dS2 

dt 



or cquivalently the corresponding evolution equations 
for the stereographic variables Wi(i) and uj2{t). Here 
Hieff and H2eff are of the same form as -ffe// given by 
eq. (3), except that the spin variable 8 in (3) is replaced 
by 5*1 and 82 in ii^ie// and if2e//j respectively, with all 
the other terms remaining the same as in eq. (3). Thus 
both the STNOs are under the influence of the common 
applied field given by eq. (6). We numerically integrate 
the above set of equations for a range of parameters in 
the /ij„ Vs a plane with random initial conditions for 
fixed values of oscillating magnetic field strength /lac and 
frequency Q, indicating regions of synchronization of pe- 
riodic oscillations (limit cycles), both in-phase (red) and 
anti- phase (black), as well as desynchronization regions 
(which includes both periodic and chaotic regimes) as 
shown in fig. [3l In fig. |3] the strong mixing of in-phase 
and anti-phase synchronizations in the (a — hA„) plane is 
essentially due to the choice of random initial conditions 
in our numerical analysis. 




hdc (O^) 



Fig. 3: The in-pliase(red) and anti-phase(black) synchroniza- 
tion regions of two similar STNOs having same anisotropy field 
K = 45.0 Oe placed in an oscillating external magnetic field of 
strength hac ~ 10 Oe and all the other parameters same as in 
fig-Ill White regions correspond to desynchronization 



Specifically in fig. 21 wc describe the synchronization 
dynamics of the spin vectors 81 and 6*2 for two identi- 
cal STNOs placed in a common external magnetic field. 
Fig. Ul^a) depicts the time series plots of signals of the 
z-components of both the spin vectors, while fig. Hl^b) 
presents the phase space structure {8l — 8l) of these com- 
ponents. Fig. m^c) displays the orbits of the two spin vec- 
tors themselves. We have chosen the system parameters 
specifically with the choice of anisotropic field strength 
K = 45 Oe. density of the spin current a = 220 Oe, dc 
component of the magnetic field, h^c ~ 130 Oe, with an 
oscillating magnetic field of strength hac — 10 Oe and 
frequency Vl = 15 ns~^ . The complete coincidence of 
the time series plots and orbits of the spin vectors along 
with the 45 ° diagonal line in the 8\ — 8^ phase space 
clearly confirms the complete in-phase synchronization of 
the magnetization vectors of the two identical STNOs. 

For a slightly different choice of parameters namely 
hdc = 106 Oe and current density a = 250 Oe with all 
the other parameters remaining the same as in fig. |4l we 
obtain an anti-phase synchronization of the two spin vec- 
tors, as depicted in fig. [51 where the y and z components 
are in anti-phase state while the x-component is in an 
in-phase state (which is allowable due to the constraint, 
8^ = 1). Note that since the direction of the spin cur- 
rent is along the y-direction (fig. ^ and the polarization 
of the spins is along the x-direction, the anti-phase syn- 
chronization along the y and z directions does not alter 
the resistance to the spin current. So both in-phase and 
suitable anti-phase synchronizations are equally desirable 
dynamical states. 
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Fig. 4: Synchronization dynamics for an array of two identical 
STNOs having same anisotropic field k = 45.0 Oe placed in 
the oscillating external magnetic field of strength hac = 10 Oe 
of frequency fl = 15 ns^^ , hdc ~ 130 Oe and a — 300 Oe. 
(a) time series plots of S], and S? (b)phase space plot of z 
component (Sz'^) of the magnetization vectors, (c) the orbits 
of the magnetization vectors S^{t) and S'^{t), showing com- 
plete in-phase synchronization. Here the open and closed circle 
points are the values of S^ and S^ at a particular instant of 
time(t ns). 



Of course real STNO devices cannot be exactly identi- 
cal and the synchronization has be to robust against mis- 
matches in system parameters to a reasonable extent. In 
order to check the robustness of the above synchroniza- 
tion of two STNOs we consider the dynamics of eq. ([7]) 
with different anisotropic field strengths, ki = 45.0 Oe 
and K2 = 46.0 Oe for the systems 1 and 2 repectively, 
with other parameters chosen as in fig. 21 The results 
are depicted in fig. [5] Again we find that the oscillators 
are in complete in-phase synchronization. In fact we have 
checked the existence of robust synchronization up to a 
choice of k2 = 55.0 Oe. 

Next we have extended our analysis to an array with 
larger number of oscillators. We have verified the above 
synchronization features for N = 4,6,10 and 100 oscilla- 



Fig. 5: Anti-phase synchronization of a system of two STNOs 
with the choice of parameters hdc ~ 106 Oe and a — 250 Oe 
and all the other parameters same as fig. [l] (a)time series plot 
of 2-components [Sl'^), (b) {Sl — Sl) phase space structure. 
(c) the spin orbits, where S^ (open circle) and S^ (closed circle) 
at a particular instant of time (t=301.5 ns) is indicated. 



tors. For N=4 in fig. [7] we have depicted the synchro- 
nization dynamics for the choice of anisotropy parameters 
Ki = 45.0, K2 = 46.0, K3 = 47.0 and K4 = 48.0 Oe with 
the remaining parameters chosen as in fig. 2] Here fig. [7^ 
presents the time series of the z-componcnts of all the four 
spins, while in fig. [7]3 the diagonal line represents the col- 
lective phase space structure {S\ — Sl~'^). In fig. [Th we 
have drawn the orbits of all the four spin vectors to demon- 
strate the in-phase synchronization behaviour in the case 
of four STNOs also. 

Finally in fig. [5] we have presented the spin dynamics 
of a system of 100 STNOs in the presence of a dc spin 
current and combined dc and ac common external ap- 
plied magnetic field of strength and frequency as in fig. |4] 



but with anisotropy parameters Ki 



1,2, ....,100, dis- 



tributed randomly between 45 — 55 Oe. Again one can 
identify a satisfactory almost complete in-phase synchro- 
nization. We have verified the synchronization dynamics 
even for a random choice of k,,'s between 45 and 80 Oe 
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Synchronization of STNOs in periodic applied magnetic field 




300 



301 f(ns) 302 



303 



-0.1 



-0.2 



-0.3 



-0.3 



-0.2 sJ -0.1 





Fig. 6: Synchronization of an array of two different STNOs 
witli anisotropy fields m — 45.0 Oe and K2 ~ 46.0 Oe placed 
in the oscillating external magnetic field, hdc = 130 Oe and 
a = 300 Oe with all the other parameters remaining same as 
fig. II 



as well, after which desynchronization sets in a slow man- 
ner (that is, a certain number of oscillators are synchro- 
nized, while the others are desynchronized) . We have also 
checked the existence of anti-phase synchronization for 
suitable choices of parameters and a combined in-phase- 
anti-phase synchronization for yet another set of parame- 
ters. The full details will be presented elsewhere. 

SUMMARY AND CONCLUSION. - In this pa- 
per we have shown the exciting possibility of synchroniza- 
tion of STNOs in the presence of a periodically varying ap- 
plied driven magnetic field in addition to constant dc mag- 
netic field and dc spin current but without any coupling, so 
that the applied magnetic field acts as a medium to induce 
synchronization of periodic oscillations. This fact coupled 
with the already proven possibility of synchronization due 
to coupling through the spin current of STNOs with or 
without delay or due to the presence of common noise, 
or mass synchronization of groups of coupled Kuramoto 
type oscillators [22], can help to make further progress 
in the experimental realization towards practical usage of 



Fig. 7: Synchronization dynamics of an array of four different 
STNOs having anisotropy fields ki — 45, K2 — 46, K3 = 47 
and K4 = 48 Oe, hdc = 130 Oe and a — 300 Oe with other 
parameters chosen as in fig. 3] 



STNOs. Further, we also note that the recent experimen- 
tal study on the action of applied magnetic field on a single 
nanowire |23] also can help the possibility in realization of 
such a mass synchronization in STNOs. 
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